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Abstract: The Witten index counts the difference in the number of bosonic and 
fermionic states of a quantum mechanical system. The Schur index, which can 
be defined for theories with at least Af = 2 supersymmetry in four dimensions is 
a particular refinement of the index, dependent on one parameter q serving as the 
fugacity for a particular set of charges which commute with the hamiltonian and some 
supersymmetry generators. This index has a known expression for all Lagrangian 
and some non-Lagrangian theories as a finite dimensional integral or a complicated 
infinite sum. In the case of Af = 4 SYM with gauge group U(N ) we rewrite this 
as the partition function of a gas of N non interacting and translationally invariant 
fermions on a circle. This allows us to perform the integrals and write down explicit 
expressions for fixed N as well as the exact all orders large N expansion. 
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The superconformal index is a beautiful generalization of the Witten index [1] 
to supersymmetric field theories. For four dimensional supersymmetric theories on 
S 3 x R it is possible to define an index which depends on three fugacities, coupling to 
conserved charges which commute with the supersymmetry generator used to define 
the index. 1 Explicit expressions in terms of integrals of elliptic gamma functions 
exist for all Lagrangian theories [2-4], 

The crucial property of an index is that it does not depend on continuous moduli 
of the theory. In the case of conformal Af = 2 theories there are many such moduli, 
the gauge couplings. These theories have a rich structure of conjectured S-dualities; 
seemingly different theories which are in fact the same under the appropriate mapping 
of couplings. Since the index should be the same for dual theories, it is possible to find 
expressions for the index of strongly interacting theories without a known Lagrangian 
(for a review, see [5]). It was also shown that the expressions for the index satisfy the 
axioms of a topological field theory in 2d living on an auxiliary Riemann surface [6]. 

A particularly simple version of the index of Af = 2 theories which depends only 
on a single fugacity is known as the Schur index [7,8]. The contributions of both the 
vector and hyper multiplets can be expressed in terms of g-theta functions, rather 
than elliptic gamma functions [9]. The topological field theory is then the zero area 
limit of g-deformed two dimensional Yang-Mills and the form of the index of the 
strongly interacting theories are given by explicit infinite sums [7]. 

The purpose of this note is to write explicit expressions in terms of elementary 
functions for the Schur index of the simplest interacting field theory, namely Af = 4 
SYM. 2 

The integral representation of the Schur index of Af = 4 SYM with gauge group 
U(N) is 3 

f> n f<j 

This expression is written in different notations than normally found in the literature; 
we use standard Jacobi theta functions and the Dedekind eta function rather than 
g-theta functions and g-Pochhammer symbols. Due to that, we also choose our g 
to be the square root of that used more often. We also rescaled a. Some simple 
definitions and relations of theta functions are given in Appendix [A]. We use the 
shorthand '&%{£) = ^(z, q ) and q = e* 7rr . 

The product in ([[]) can be rewritten using an elliptic determinant identity. For 

4 In addition one can define fugacities for flavor symmetries. 

2 Other limits of fugacities of the index of Af = 4 SYM were studied in [10]. 

3 For SU(N) gauge group the prefactor should be rescaled by g 1 / 4 r7 2 (r/2)?y _4 (r). In addition a 
delta function should be introduced so that sum a* vanishes modulo ir. The final answer is the same 
as for U(N ), up to the aforementioned prefactor. We thank Arash Arabi Ardelrali for pointing out 

an omission in this factor. 
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arbitrary cp and a', one has [11-14] 
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Here we used the notation = $*(0, q). Setting a[ = cni simplifies the prefactor and 
we can also use the quasi-periodicity of $3 to write 


$ 3 (tttN/ 2) =q n2/4 V 3 A n , 


Ajy — 


N even, 


JVodd. 
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Now we note that the ratio of Jacobi theta functions appearing in the determinant 
is in fact closely related to the elliptic function cn with modulus k = d\l"d\ 


= ( 2) 
tf 4 (z) ’ 

Finally using the relation = 2rf{r) we can write the index © as 4 
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Equation (Q) has the form of the partition function of N free fermions on a circle 
(see e.g. [15]). The Fermi gas partition function is completely determined by the 
spectral traces 


Z(_ = Tr(//) = / da 1 ... dai p (aq, a 2 ) ■ ■ ■ P ( a e, a i) , 

Jo 

q 2 1 i(2p-l)(a-a') 

>>(“.“') = 5 Jcn((a-Oi')^) = -J] 


( 6 ) 


q p 2 4 - q 


-( p - 2 ) 


where we used the Fourier expansion of the cn function. 

Performing the integrals over a identifies the Fourier coefficients in the expansion 
of the different p’s giving an exceedingly simple result 

t 




rP-§ + q -(P~h) 


(7) 


As we show in Appendix |B|. these infinite sums can be expressed as polynomials 
of complete elliptic integrals K and E. One can then use the combinatorics of the 
conjugacy classes of the symmetric group Sjy with mg cycles of length t to write 
Z(N) (|) as 

mA 

M e 1 


z w=£'ii 


4 The factor of $|/27r is included with the cn functions to simplify their Fourier expansion below. 
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where the prime denotes a sum over sets that satisfy Y2e & m e — AT. Plugging ( |B..5| ) 
into (||) and including the normalization in (j5j), we find for N — 1,2, 3,4 

„-i/4 cr l 

X(l) = q - - VkK , 1(2) = ^~ 2 K(K - E ), 

7r Z7T z 

— 9/4 

1(3) = -VkK(l2K(K — E) — 4(1 + k 2 )I< 2 + vr 2 ) , (9) 

247 T 6 

1( 4) = J-2 -K(3K(K - E) 2 - 2 k 2 K 3 + n 2 (K - E)) . 

It is easy to generate explicit expressions for larger N. The powers of q in these 
expressions as well as y/k for odd N come from the normalization prefactor in (|5|). 
We see that by removing this power of q, the resulting expressions are polynomials of 
complete elliptic integrals and the elliptic modulus and have nice modular properties. 
We comment on that further below. 

An alternative approach is to consider the grand canonical partition function 
with fugacity n 

OO 

2 («) = 1 +£>( JV > t "- ( 10 ) 
N=1 

This is a Fredholm determinant of a very simple form 


S(/c) = exp 



=n 
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This product turns out to be expressible in terms of theta functions (see Appendix |A|) 
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Recall that the index has the extra factor of 1 /Ajv (0) compared with the free 
fermion partition function (^]). This factor distinguishes even and odd N, so it is 
natural to split S into its even and odd parts. In fact, since $ 3 (z) is even and d 2 (^) 
is odd under z —» 7 r — z, the two terms in the last line of (|l 2 j ) give the decomposition 
into 

S ± = -(s( K )±S(- K )) (13) 

It is then possible to define the grand index S, where the EL part is divided by the 
value of A n for odd N 


E(k) = 1 + ^l( N )q N2/A n N 
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By expanding the last expression in powers of k, it is very easy to obtain explicit 
formulae for X(iV) for finite N. Using relations between derivatives of theta functions 
and elliptic integrals, one recovers 

We can also obtain a closed formula for the coefficients of the large N expansion 
of 1{N) by expanding ([14]) at large k = e ;i and using the integral transform 


Z(JV) = 


j-N 2 /4 
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dn H(e M )e 
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Given the large k expansion of arccos 


arccos — = i arccosh — = i/j, + i log 
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where the exponent is normal ordered. Expanding in powers of e 2/i gives 
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The index (|U|) can now be easily evaluate through integration by parts. For any 
integer l > 0, we have 

\l riir 


— [ du e-^ N+2n) d l J 3 (in) = ( N + 2n ) [ dne-^ N+2n ^ 3 (in) 
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and similarly for $ 2 (*/•*) • Using the Fourier expansion 


I dfie ^ N+2n) +t? 2 (*//)) = q (N+ * n) , 


we finally obtain 




n=0 


AT + n 
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( 20 ) 


( 21 ) 


This sum is convergent for \q\ < 1 and it is easy to check that for N = 1, 2, 3,4 it is 
exactly the g-expansion of the finite N expressions in ([|). 

Equations (|1]), (|T4|), (|12j) and ([|) (with the algorithm in Appendix jB{) are the 
main results of this paper. 

The large N expansion (^lj) should have an interpretation in string theory, but it 
is a long enduring puzzle to get a match. The leading contribution to the Schur index 
with n — 0 is just l/'d 4 , which is N independent, unlike the supergravity action, which 
should scale with N 2 . This was already observed for the large N limit of the most 
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general index of AT = 4 SYM in [3]. Indeed, the correct quantity to match should be 
the partition function on S' 3 x S 1 , which is related to our definition of the index by the 
Casimir factor, which does scale as N 2 [16,17]. It still does not match the gravity 
calculation, possibly due to missing counterterms of supersymmetric holographic 
renormalization. 

Still one could hope that the terms with n > 1 correspond to some instantons 
in string theory. Identifying logl /q with the radius of the compact S 1 (relative to 
the radius of S 3 ), it seems to match the action of n D3-branes plus n 2 excitations 
between the branes. Natural candidates would be 1/8 BPS giant gravitons [18-20], 
which should play a similar role to the string and membrane instantons in the case 


of ABJM [21,22], 

The grand canonical partition function 5 (|l^) and the grand index S (14) bear 
a very strong resemblance to the grand canonical partition function of Af = 8 super- 
conformal Chern-Simons-matter theories, which was expressed in [23] also in terms 
of a Jacobi theta function. In the 3-dimensional case, the classical limit of the grand 
potential for a wide class of theories contains a cubic term in p which leads to the 
universal Airy function behavior for the partition function [24], The analog in 4 
dimensions is a quadratic term in p (related the the structure of the Fermi surface 
on the circle). The analog of the Airy function is then a q kN2 scaling of the partition 
function of the fermions, here with k = 1/4. But this behavior is not evident in the 
final answer, due to the overall prefactor in the definition ([!]), ( 0 ). If the p 2 scaling 
is universal, it would be natural to try to reproduce it from a quantum supergravity 
calculation, along the lines of [25], though an understanding of classical supergravity 
is required first. 

Comparing to the Fermi-gas formulation of 3d theories [24], the expressions we 
find are even simpler, since the density p ( 0 ) is a function of difference form, which 
means that the Hamiltonian H = — log p is a function only of the momentum. So the 
free fermions have a complicated kinetic term, but no potential and no interactions. 
Indeed the exact spectrum is pairwise degenerate with 


E P = log (q p 5 +q = - (p -^ logq - ^ i), P > 1 • (22) 

Identifying h = —inr, this is the spectrum of the harmonic oscillator with exponen¬ 
tially small corrections at large h. In the 3d case there is a connection between the 
grand canonical partition functions and the spectrum of Schrodinger equations on 
the line related to topological string theory [26,27]. One can easily enrich our model 
such that it would calculate non-trivial spectra on the circle. A natural modification 
is to include extra fundamental hypermultiplets. While the resulting matrix model 
looks interesting, the 4d field theory is not asymptotically free, so it is not obvious 
at all what this would be capturing. 
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It was advocated in [9] that the index would have nice modular properties once 
rescaled by g c , where c is the central charge. This rescaling is exactly that mentioned 
after (|9j) and is also the one appearing in our definition of the grand index 5 0)- 
Indeed for finite N we find elliptic integrals and at fixed /x we get theta functions 
of norne q. all of which have nice modular properties. It was also proposed in [9] to 
modify the index by changing some charge assignments, which amounts in our case 
to introducing a mass m — 7r/2 for the adjoint hypermultiplet. It is easy to repeat 
our calculation for this particular mass and the result for the grand index is equal 
to our S_|_(ze M , — q ) (jT3j). The modified index therefore vanishes for all odd N, and is 
essentially the same as the usual index for even N. 

It would also be interesting to reproduce our expression from the topological 
held theory formulation as a sum over U(N) representations [6,8]. Given the explicit 
relation in this case to g-deformed 2d Yang-Mills [7], it should also be possible to 
reproduce our results from an analysis similar to [28]. The Schur index is also con¬ 
jectured to be equal to the torus partition function of certain 2d chiral algebras [29]. 
For AT = 4 SYM the conjectured chiral theory has super W-algebra symmetry and 
our expressions should give the full partition function of this theory. In addition 
there is a relation between the index and integrable models [30], which should also 
yield the same result. 

One can also apply such ideas to other theories: Does the grand index take a 
simple form for all class-iS theories? The most obvious generalization is necklace 
quivers, where a lot of the structure persists. It is also possible to include fugacities 
for the flavor charges, which in the case of AT = 4 is analogous to considering Af = 2*. 
Other generalizations are the inclusion of line operators, for which the matrix models 
were computed in the Schur limit in [31]. It would also be interesting to see whether 
it is possible to generalize our analysis beyond the Schur limit and/or to the lens 
space index [32], We hope to report on some of these issues in the near future. 
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A. Theta functions 


In the literature on the superconformal index, the expressions are mostly written 
in terms of g-theta functions and q -Pochhammer symbols. Those are related to the 
Jacobi theta functions and Dedekind eta function by 

(«=*: = n (1 - «*) = - (A.1) 

where the (quasi) period r is related to the nome q by q = e™ r . 

The Jacobi theta function i9 3 (z, q ) is given by the series and product representa¬ 
tions 


oo oo 

Uz,q)= V ?”V i ”' = n(l-?“)(l + 2? 2t - 1 co S (2z) + I! "- 2 ), (A.2) 

n =—oo k =1 

in terms of which the auxiliary theta functions are 

tfi (z, q) = iq l/i e~ lz d 3 (z - r - ±7T, q) 

02 (z, q) = q l,A e~ lz d 3 {z - tt, q) (A.3) 

04 (z,q) = 0 3 (z - T,q ). 

We also need Watson’s identity (eq. 20.7(v) in [33]) 

<l 3 (z,q m ) 03 (w, 9 1/2 ) = Uz + w, <?)0 3 (z — w,q)+ ^(z + W, g)0 2 (^ - w, q) (A.4) 


B. Finite N expressions 


The quantities Z\ (jTj) can be evaluate by the procedure given in [34] which we briefly 
summarise (correcting some small typos). For even 1 = 2s one has 

Zls = 2 ^ ^ a m (s)P2s-2m-l ) (B.l) 

' ' m =1 


where a m (s) and D s are generated by 


s-i s -1 00 / 4/2\ 

£>»,(»)<"■ = ip 

m— 0 j= 1 s =0 ' 


EK 

2tt 2 


K 2 
2 t r 2 


(1 — k 2 ) nd 2 



J 


(B-2) 

where AT = K{k 2 ) and E = E(k 2 ) are complete elliptic integrals of the hrst and 
second kind respectively, nd is an elliptic function (likewise cd below) and the elliptic 
modulus is k — ‘^/d 2 . 
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Similarly for odd i = 2s + 1 


Z 2 s+i — 


- 1 )‘ 


2 2s_1 (2s)! 

v J m =0 


^ ^ (^) J c * 


2s—2m 1 


where 


£<u»)« m =n i 1 - ( 2 j - i > 2 d - 

m= 0 j=l 

We find for N = 1, 2, 3,4 
kK 

Zi = —, 

7 r 


Z, = — 


Zi (1 - fc 2 )fcW 3 


27T 3 


(2s)! 2vr 


- (1 - k 2 )K 2 


Z, = 


n- 


Z 2 (1 - k 2 )k 2 K A 


Z 4 = — 
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3t r 4 


(B.3) 


y. , (~t 2 )’ Kk , f 2Kt \ 
/ J2s~— rr _ cd 

s=0 


7T / 


(B.4) 


(B-5) 


Plugging these expressions into (^) and including the normalization in (|5|) gives (p). 
This algorithm can be easily pushed to larger values of N and is easy to implement 
in Mathematica. 
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